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Abstract 

The deuteron disintegration at high energies and large angles in the 
d{p, 2p)n reaction, is calculated in kinematical conditions where the dom- 
inant contributions are due to soft rescatterings of the initial and final 
nucleons, which accompany the hard pp reaction. The eikonal approxi- 
mation, which accounts for relativistic kinematics as dictated by Feynman 
diagrams, reveals the important role played by the initial and final state 
interactions in the angular and momentum dependences of the differential 
cross section. Based on these results, we propose a new and effective test, 
at moderate energies, of the physics relevant for the color transparency 
phenomenon in hadron-initiated exclusive hard processes. 



I. INTRODUCTION 

We have demonstrated recently [jl| that it is possible to regulate the relative contri- 
butions of different reaction mechanisms in the exclusive high energy hard pd ppn 
reaction, by selecting specific kinematical conditions. In that paper we discussed this 
reaction within the kinematics that is preferential to the study of short range nucleon 
correlations in the deuteron, the role of relativistic effects in the deuteron structure and 
the effects of nucleon binding. The aim of the present paper is to analyze the reaction in 
the kinematics favorable for the search of the Color Transparency (CT) phenomenon. 
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It has been shown in refs. |T|J^ that, under certain kinematical conditions, initial 
and final state interactions can dominate the cross sections. Those same kinematics 
can be advantageous for studying CT by the vanishing initial (ISI) and final (FSI) state 
interactions of hadronic products in the hard exclusive pN scattering. (For details on CT 
physics, see Ref. and references therein). As a first step for the identification of CT 
it is necessary to develop a baseline model which accounts for the initial and final state 
interactions within the conventional eikonal approximation. Thus, in this paper, we first 
restrict our considerations by calculating the elastic rescatterings of nucleons only. Next 
we will include the inelastic production in the intermediate states in terms of CT. 

In Ref. the high momentum transfer exclusive deuteron electro-disintegration reac- 
tion has been calculated within the eikonal approximation. The calculation of the Feyn- 
man diagrams, relevant for the propagation of knocked out energetic nucleons through 
the nuclear medium, revealed certain effects related to the finite longitudinal distances, 
which are beyond the formulae of conventional Glauber approximation One of the 
aims of this work is to extend the formalism of Ref. to account for similar effects in 
the d{p, 2p)n reaction. 

Another motivation of our study is to understand the unexpected energy dependence 
of the semi-inclusive A{p,p)X reaction. In the energy range up to the 10 GeV, Carroll et 
al. Ref. P found that, in semi-inclusive A{p, 2p)X reactions, the absorption of protons 
was reduced compared to the value predicted in the Glauber approximation. At incident 
proton energies larger than 10 GeV an increase in the absorption was observed. That 
result led to a number of theoretical suggestions, such as the possible importance of the 
interference between interactions of point like- and large size configurations ||^, the longi- 
tudinal momentum dependence of nuclear transparency at pre-asymptotic energies [pO| , 



the possible existence of a threshold for charmed particle production pj], etc. Explor- 
ing the deuteron in the above energy range and restricting the kinematics to where the 
deuteron wave function is well known, may help to elucidate the physics of hadron prop- 
agation through the nuclear medium and the energy dependence of the semi-inclusive 
reaction. 

As a definition for hard pp collision in the exclusive d{p, 2p)n reaction we choose the 
conditions of large center of mass angle (> 60^) pp scatterings with |m| > 2 {GeV/c)"^. 
We restrict the calculations to momenta of the incoming proton pi > Q GeV/ c, where a 
reliable parameterization for hard pp scattering exists (see e.g. |l|). 

We deduce eikonal formulae of the scattering process by analyzing the Feynman di- 
agrams corresponding to the single and double rescatterings(see Fig.l). At the present 
kinematical conditions the nonrelativistic deuteron wave function has been used but the 
full relativistic kinematics of the interaction is taken into account by the Feynman dia- 
gram approach. At intermediate energies we account for relativistic kinematics exactly, 
i.e. we keep also terms ~ whose contributions are enhanced due to the steep mo- 
mentum dependence of the deuteron wave function. We make the approximation that 
the momentum of the fast nucleon is almost conserved in elastic high energy soft (small 
angle) collisions. Thus, the hard process kinematics are not strongly infiuenced by ISI 
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and FSI (see Appendix B). 

In section 2 we outline the details of the calculation, and explain the difference between 
our results and the nonrelativistic Glauber approximation. The results can be used to 
calculate nucleon rescatterings off any nucleus. For the deuteron we calculate analytically 
the scattering amplitudes, using the known (see e.g. Refs. [0,|T^,|I^) parameterization 
of the wave function in momentum space. 

In section 3 we present the results of numerical calculations for the kinematical condi- 
tions that can be reached in the EVA experiment at BNL [^. The results reveal a strong 
impact of initial and final state interactions on the angular and momentum dependences 
of the differential cross sections. A surprising prediction is the azimuthal dependence 
of the final state interaction, which is a consequence of the interplay between hard pp 
scattering and soft pn rescatterings. We analyze also the range of applicability of the 
commonly used factorization of hard and soft scattering amplitudes. In section 3 we es- 
timate the size of expected CT effects and suggest new options for the search for various 
consequences of CT. The simultaneous observation of all effects will make it possible to 
reach unambiguous conclusions about the onset of CT. The results will be summarized 
in section 4. 

In Appendix A the analytical calculation of the rescattering amplitude for the deuteron is 
presented. In Appendix B we estimate the distortion of the kinematics of hard processes 
due to soft A^A^ rescatterings. 



II. KINEMATICS AND APPROXIMATIONS. 

In this section we calculate the cross section for the d{p, 2p)n reaction in the kinematics 
where the deuteron is at rest and the pp scattering is "hard". The kinematics of the 
reaction is determined by the four-momenta pd = {Ed,Pd), Pi = {Ei,pi), p^ = (E^yPs), 
P4 = (£^4,p4), and Ps = {Es,ps) of the deuteron, incoming and scattered protons and 
spectator neutron, respectively (c.f. fig. la). We define also the variable as = ^^^12 ' 
which is convenient for the description of high energy processes. The quantity ^ can 
be viewed simply as the fraction of the deuteron momentum carried by the spectator 
neutron, in the infinite momentum reference frame ( in which the deuteron is fast). The 
indices "t" and "z" denote the transverse and longitudinal directions with respect to 
the incoming proton momentum p^. The light cone momentum of the target proton is 
a = 2 — as- We will use also the spectator azimuthal angle (ps, which is the angle between 
the (pi,P3) and {pi,Ps) planes. 

We choose \t\ = \{pi — ps)^] ~ izi^n- and \u\ = \{pi — Pa^I > 2 GeV"^ in order 
to fulfill the requirements of hard pp scattering at large angles. With these kinematics, 
in the rest frame of the nucleus, the transverse components of the proton momenta are 
small compared to the longitudinal components i.e.: 

'^3^^2 _ ^2 a< 1; (^)2 _ ^2 _ 2— (2 - a) < 1 (1) 



P^z Pi Piz Pi 
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where 6^ and 6^ are the polar angles between the scattered/ knocked-out protons and 
the incoming proton (pi). 

We restrict the kinematical region to relatively small spectator neutron momenta i.e. 
as = 2 — a ~ 1 and pst < 300 — 400 MeV/c. It was shown in ref.( [jl|) that vacuum 
(Z)-type diagrams in the d{p, 2p)n reaction can be neglected under those conditions. 
By restricting the calculation to this range of Fermi momenta we can also neglect the 
contribution from possible non-nucleonic components in the deuteron wave function. As 
a result it is safe to use a conventional non-relativistic deuteron wave function (see e.g. 



Refs.f [12,13 



Fig.l depicts all the relevant Feynman diagrams for the impulse approximation and 
for the initial and final state interactions, when the soft rescatterings are calculated within 
the elastic eikonal approximation. The solid circles in the diagrams represent the hard pp 
scattering amplitude and the broken lines represent the amplitude of small angle elastic 
NN scattering. Note that, in eikonal approximation, there can be no contribution from 
diagrams where a fast proton experiences soft rescatterings off both target proton and 
slow spectator neutron. This is due to the finite longitudinal distance between the target 
proton and neutron and the fact that the sequence of the projectile proton colliding 
off the target proton and neutron, followed by a hard pp scattering is geometrically 
impossible. Similarly, it is easy to demonstrate that only diagrams of Fig.l survive 
within a eikonal approximation. A similar cancelation of diagrams is well known in 
nonrelativistic quantum mechanics where infinite series over exchange by potential are 
reduced, in the high energy limit, to the diagrams expressed through the full amplitudes 



of NN scattering in 



A. The impulse approximation (lA) 

We can classify the Feynman diagrams in Fig.l by the number of soft initial and final 
state interactions. The simplest amplitude is the impulse approximation amplitude of 
Fig. la. that can be written as: 

F, = {27T)h{p:,p,)A'^^^%,t), (2) 

where A'^"'^'^ is the amplitude of the hard pp scattering and ip is the nonrelativistic 
deuteron wave function normalized as / ip'^{p)(Pp = 1. We are interested in the cross 
section on a deuteron normalized to the pp cross section at the same s and t, so that the 
hard pp scattering amplitude is cancelled out. 



B. The single rescattering amplitude 

In the eikonal approximation the first order rescattering is described by the diagrams 
of Fig.lb,c,d. The single rescattering amplitude ( Fig. lb) can be parameterized as 
follows: 
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,2 _ ^2 I ^Vl r^/2 _ ^2 I ,Vl [/'^ . _L ^ _ ^/ ^2 _ ^2 I ,Vl ,Y0-Tr^4 ' 



2m [pI — rn? + ie] [p'^ — m? + ie] [{p^ + Ps— p's)"^ — m? + ie] i(27r) 



(3) 



where, in the intermediate state, p'^ is the momentum of the spectator neutron, s' = 
{P3 + P4, + Ps — p's)^ , ^{Pd,Ps) is the invariant vertex of the transitions d ^ pn into two off- 
mass shell nucleons and /^^^ is the spin dependent amplitude of pn soft scattering. All spin 
dependences of the target nucleons is included in the vertex factor. We factorize the hard 
scattering amplitude A^"^'^ out of the integral because the transferred momenta in soft 
rescattering are negligible compared to the transferred momentum of the hard amplitude 
(The validity of this approximation is discussed in Appendix B). It is reasonable to 



neglect also the dependence of fpn on p'^Q (see c.f. [|15|) because, at high energies, the 
NN interaction depends only weakly on the collision energy. Using a non-relativistic 
description of the Fermi motion in the deuteron allows us to evaluate the loop integral 
by taking the residue over the spectator nucleon energy in the intermediate state i.e. we 
can replace [p']' — rn? + ie\~^d^p'g by —i{2n)/2E'^ ^ — 1(271) /2m. This is possible because 
in this case it is the only pole in the lower part of the p'^^ complex plane. The calculation 
of the residue in p'^^ fixes the time ordering from left to right in diagram Fig. lb. We 
introduce the nonrelativistic deuteron wave function as: ipiPd — p\) = — -7 — 5 — 

(with / |?/'(fc)p(i^fc = 1) (see e.g. [llSHr^)- In the laboratory system the amplitude can 
be rewritten as: 



{^T^) '^_j^hard,^, ^^ I „u P^'iPi.P's.Ps) d^P's 



2m J ^ [{p4 + Ps — Pg)"^ — m"^ + ie] {2tt) 
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Wi,hard(^ ,^ f /^ f^Psf-p'st) d^P's 

2 ^ ^^''V^^ ^^^2mp4.[j>L-J>.. + A4 + .e](27r)3' 



where 



A4 = {Es - m)— - {pst - p'st) — - (5) 

P4z Piz 

In the last part of eq.(^) we used energy-momentum conservation to express the 
propagator of a proton of momentum p^ as: 



(P4 +P3- Psf -m? + ie 

'^.pAz 



, fj^ . , .Pit , (Ps-P's)^ 

Psz - Psz + {Es-m) [pst - pj — + 



P4z P4z "^PAz 

"^Pizlp'sz ~ Psz + ^4]- (6) 



We use the fact that the energy transferred in the pn rescattering is negligible compared 

2p4z 



to the total energy of the scattered particles so that we may neglect the term: ^^-^ 



with respect to the other two contributions to A4. 

We keep the term (Es — fn)j^ because it does not vanish with the increase of the 
projectile energy at fixed spectator nucleon momentum. We also keep the term [pst — 
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p'st)^ which decreases only as according to eq.(|I]). If the projectile energy is not 
extremely large, the contribution of this term in some kinematical conditions can be 
comparable in size to the first term of eq.(|^). 

The fact that the soft A^A^ scattering amplitude depends only weakly on the initial 
energy helps to simplify eq.(^). It is convenient to redefine the soft scattering amplitude: 
2^4' m'^"'' ~ /^"(Ps ~ Ps)) where /p" is now the scattering amplitude normalized by 
the optical theorem /^"(fc = 0) = icrfot- We may neglect the longitudinal momentum 
transfer in f^^{p's — Ps) and assume that the transferred momenta are almost transverse: 
i.e. Ps — p's ~ Vst — v'st = ^t- Approximating the transferred longitudinal momenta, 
as kz = Psz — v'sz — ^4 using the fact that in the rescattering integral the average 
transferred momenta in the pn scattering are ~ pt, we see that the condition A4 <^ \pt\ (or 
\ps\ ^ 2m) allows us to neglect the longitudinal momentum transfer in the soft scattering 
amplitude. Therefore, under our kinematical conditions, eq.(§) can be rewritten as: 



2 J [p'^^ - Psz + A4 + ie\ (27r)'^ 

The integration in eq.(|^) can be performed in coordinate space. Writing the Fourier 
transform of the deuteron wave function as: 



m = 7— T / d'r<P{r)e-'^^ (8) 
(27r)2 

and using the coordinate space representation of the nucleon propagator: 

^ -i I dz°e(z°)e'(P'-~P-+^^)"° (9) 



[P'sz - Psz + A4 + ie] 
we obtain the formula for the rescattering amplitude: 

= -v4'^'^^'^(s,t) J rfV0(r)0(-2)rP"(A4, -z, -6)e^f«'' 
where r = — r„, kt = Pst — p'st- We define a generalized profile function F: 

FP'^(A4, h) = ie-^^*^ j r{kt)e-^^^'0^^. (11) 

Eq. ( pTD reduces to a Glauber- type approximation in the limit of zero longitudinal momen- 
tum transfer A4. The dependence of the profile function on the longitudinal momentum 
transfer originates from two sources. According to eqs.(^,|^) the term [Eg — m)^ ac- 
counts for the relativistic kinematics which follows from the evaluation of the propagator 
of a proton with momentum ^4. A similar factor was found also for the c?(e, e'p)n reaction 
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in ref. The second term (pst — p'st)^ — is due to a slight misahgnment of the 
projectile momentum with the directions of the outgoing protons. Even though, accord- 
ing to eq.(|l]), this factor is small it has a large effect on the cross section due to the steep 
momentum dependence of the deuteron wave function (see discussion in section 3). The 
same modified profile function, eq.([Tl|), is valid for the single rescattering amplitudes off 
any nucleus A [0. The unique feature for the deuteron is the possibility to calculate 
the rescattering amplitudes using an analytical parameterization of the deuteron wave 
function (c.f. |2|,|T3l). 

Substituting the deuteron wave function in eq.(|^) by the analytic form of eqs.( |Al|) , 
and (|A3D, one can perform the integration over p'^^ by deforming the contour of 



integration in the upper complex half-plane of p'^^. We can then evaluate the integral by 
taking residues over corresponding poles in the deuteron wave function. The answer (see 
appendix A for details) for the single rescattering amplitude Fb (eq.(^) is: 



Ai ' ' '7 (27r 



^ (^''^) / TTT^.n^,) iriPs) - #'^(Ps)) , (12) 



where Ps = Ps{Psz — ^4, Pst — kt), is deuteron wave function defined in eq. (|Al|) and 
is defined in eqs.([A7D and ([A^). 

We can now calculate also the amplitude corresponding to the diagram of Fig.lc, 
replacing p^ — > Ps- Thus 



A^'^'^s.t) / rfV0(r)^(-2)rP"(A3, -6)e*P^ 



(2^)lA'^-'^(,,t) / ^,r{kt) iriPs) - li^'^iPs)) , (13) 



Ai ' ' (27r 

where ps = psipsz — ^3, Pst — kt) and A3 is described by eq.(||) with the substitution 
P4 Ps- 

The initial state interaction represented by diagram Figl.d. is described by an equa- 
tion similar to eq.(|^: 

F_ (27r)^ ,^ /• /^ /"""(Pst-P'st) d^P's ...x 

where Ai is defined by eq.(|]) with the substitution p^ ^ pi. In contrast to eq.(^ the 
singularity of the propagator in eq.(|T^) is now located in the upper part of the complex 
p'g^ plane. It is easy to show that, in coordinate space, this difference corresponds to the 
reverse sign of the argument of 6 function in eq.(|^): 0{—z) — > 9{z). In momentum space, 
the contour of the p'^^ integration can be deformed in the lower complex semi-plane, where 
the only pole is due to singularities in the deuteron wave function at p'^^ = —i^Jpft + 
(see Appendix A) and the result of the integration differs from eq. ([T2|) by the sign of p'^^ 
only. Thus the amplitude for diagram Fig. Id is given by: 
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Fa= J d^r(j){r)e{z)TP''{Ai,-z,-b)e'P' 



Ai ' ' '7 (27r 

— * 

where ps = Ps{Psz - ^i,Pst - h)- 



C. The double rescattering ampUtude 

The second order rescatterings where both protons scatter off the spectator neutron 
are represented in the diagrams of Fig. le,f,g and h. Using the same approach as for the 
single rescattering, the double rescattering amplitude of Fig.le can be written as: 



ipi-pl 



[pi'J - p^) - Ai + le] [pf] - Ps. + A3 + le] (27r)3 (27r)3 

(16) 



where p^p and p^f^ are the spectator (neutron) momenta before the first and second 
rescatterings, respectively. 

Since the deuteron wave function in eq.(0) does not dependent on p'f\ the whole 
dependence on pf) is contained in the two propagators only. But the poles of these 
propagators in the pf) complex plane are located in the lower semi-plane. Due to fast 
convergence of the integral, the contour of integration can be moved to the upper complex 
semi-plane which makes the integral equal to 0. In the coordinate space this result has 
a simple geometrical interpretation. The diagram describes an interaction of incoming 
and scattered protons with a spectator neutron which would be located simultaneously 
before and after the target proton which is not possible geometrically. Indeed using the 
coordinate representation of the deuteron wave function (eq.(||)) and the nucleon prop- 
agators (eq.(^), we can obtain the Q{z)Q{—z) term in the integrand in the coordinate 
space. The same reasoning leads also to zero for the amplitudes of diagrams Fig l.f, 
which can be obtained from diagram Fig l.e by substituting p^ p^. The amplitude in 
diagram Fig. Ig, within the approximations discussed above takes on the form: 



(27r)i f g). r"(pi? - P^t^) f^Pst - pg) d^^^ d^pf) 

^'^ ' [pii^ - pi'J + A3 + ^e] [p?J - Ps. + A4 + ^e] (27r)3 ^ > 

where A3 = {E^ - m)|^ - (pif - p«)g and A4 = {E, - E^^))^ - {p,, - pi?))£ 
are different from those for single rescattering (eq.(|^)). However, in our kinematics 
Es - Ef^ ~ ^ < m and A3, A4 can be taken to be the same as for single rescattering. 
Using eqs.(H) and eq.(|^) one can write the amplitude, in coordinate representation, as 
follows: 
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Fig) = \j rfV0(r)^(-z)r?'"(A3, -6)r^^"(A4, -z, -6)e*^^^ (18) 

where F^" is defined by eq.(|lT]). The difference between eq.(|18[) and the conventional 
nonrelativistic Glauber approximation is the presence of the terms A3, A4. The integra- 
tion in eq.(p!7|) over the longitudinal components of nucleon momenta p^^] and pf) in the 
intermediate states can be evaluated by the technique used in Sec. 2. 2 for the calculation 
of the single rescattering amplitude(see Appendix A). First, we integrate over pf) by 
closing the contour of integration at one of the poles in the complex plane (they are 
located now in different semi-planes) to obtain: 

8 [pVJ - Psz + A3 + A4 + te] (2vr)2 27r (27r)2 

The p^^J integration in eq.(0) can be performed in a way similar to the one used in 
eq.(l^) for A = A3 + A4. Therefore, we may use the result of the calculation of single 
rescattering to obtain: 

Fg = / ^^rik^'^)riki'^) mps) - ^^''{Ps)) , (20) 



16 J (27r)2 (27r) 

where ps = Ps{Psz — A4 — A3,pst — kf'^ — kf'''). The expression is symmetric under 
interchange of A4 A3. Thus, for the amplitude of diagram Figl.g we get: 



Fg (21) 



D. The cross section in DWIA 

In the approximation where the hard pp scattering amplitude is factorized (Appendix 
B) the cross section of the d{p, 2p)n reaction can be written in a form which is similar to 
the distorted wave impulse approximation (DWIA). Hence the cross section for d{p, 2p)n 
reaction can be expressed as follows (see e.g. Ref. [^): 

n{p,)b{E, - {Md - E)) (22) 



{d^p^/E^){d^pjEi) 271 2m-\pi\ dt 

where E = E^ + E^ — Ei. The function n{j)s) describes the distorted momentum depen- 
dence of a nucleon in the deuteron: 

_ \F. + F, + F, + F, + Fg + F,\- 
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The factor (27r)^ follows from the phase space factor of the spectator nucleon wave func- 
tion. In eq(^2D we neglected the virtuality of the interacting nucleon [| and used for 
the experimentally measured cross section and for which a phenomenological parameter- 
ization is given in Refs. [p^flfl. 



III. NUMERICAL ESTIMATES 

We wish to calculate the quantity T- "Transparency" which accounts for the effects 
due to soft rescatterings. Traditionally (see e.g. [^) T is defined as the ratio of the 



experimentally measured cross section to the calculated lA cross section, which according 
to our definitions is: 



^ ^DWiA ^ \Fa + F, + F, + F, + F, + Fn\ 



a 



For our numerical estimates we will use the kinematics of on-shell 90° center of mass 
angle pp scattering. This is achieved by imposing the condition 



Aw? — s 

t = (25) 



on 



s = {P3+ Pif and t = {pi- p^y-. 



A. Calculation within the elastic eikonal approximation 

We shall call "elastic eikonal", the approximation where nucleons only propagate 
through the deuteron and soft pn rescatterings are described by the diffractive f^^ scat- 
tering amplitude. Since in the relevant kinematics the amplitude is predominantly 
imaginary, the lA (eq.(0)) and double rescattering amplitudes (eq.(^)) are positive, while 
the single rescattering amplitudes (eqs(|l^), (|T3|) , (|15D ) are negative. Therefore the terms 
in eq.(^) have different signs and, depending on the kinematics, the soft rescattering 
terms can either increase or decrease the transparency T, leading to large ISI/FSI effects 
on the overall cross section. In Refs. we demonstrated that the ISI/FSI is largest 
for perpendicular kinematics, where the polar angle of the neutron-6's is almost perpen- 
dicular to the reaction axis. In Fig. 2 we show the spectator angular {6s) dependence of 
the transparency T, for different spectator momenta ps- The calculation was done by 
adopting the standard parameterization for the diffractive pn scattering amplitude: 

/^" = c^fo"(^ + «n)e-^"'=?/^ (26) 



This approximation is legitimate as long as we neglect the production of intermediate inelastic 
states. 
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In our kinematics, the parameter values 0"^! ~ 40 mb, an ~ —0.2 and 6„ ^ 8 GeV ^ were 



obtained by fitting the experimental pn scattering data |24]. The rather weak dependence 



of atot and 6„ on the nucleon momentum, for energies above 3 GeV, has been included. 

The results, presented in Fig.2, confirm the large ISI/FSI effects in perpendicular 
kinematics. They reveal also the importance of the above discussed modifications in the 
conventional Glauber approximation |Q . We find that the dominant effect is due to non- 
zero transverse momenta of the emerging fast protons (~ kf^^). This effect reduces 

strongly the contribution from the interference of two rescattering amplitudes. In Fig. 3 
we show the pt dependence of T at fixed value of = 1. It follows from this figure 
that in the range pt < 200MeV/c the ISI and FSI significantly screen the cross section. 
It is small compared to the cross section calculated within the lA approximation. The 
interference between single rescattering and the lA term dominates the overall scattering 
amplitude. With increasing > 250 — 300Me^/c the squared terms and the interferences 
between single rescattering amplitudes dominate and as a result T may exceed 1. At even 
higher pt the second order (negative) terms in eq.(P^) (e.g. interference between single 
and double rescattering terms) tend to suppress the rise of T, thus changing the slope of 
the Pt dependence. 

Since at fixed 9c.m., 6*3 and 64 depends on the incoming proton energy (see eq.(|I])), the 
formulae of eikonal approximation could lead to some energy dependence of T for fixed 
spectator momentum, even if afot, bn were energy independent. This dependence is more 
noticeable for the larger momenta of the spectator, since the contribution from ISI/FSI 
diagrams increases with increasing pt (see Fig. 4). In fact, we see in Fig.4 that T does not 
depend much on incident beam energy. 

Another nontrivial consequence of the interplay of hard pp scattering and soft pn 
rescattering is the dependence of FSI on the azimuthal angle of the spectator neutron. 
The effect is due to two important features of the reaction: in hard pp scatterings protons 
are produced at small but finite 9^ and 64^ angles (lab) and high energy soft pn rescattering 
is characterized by an average transferred momentum < k"^ >^ pf, which is perpendicular 
to the trajectory of the protons. To visualize the origin of the azimuthal dependence let 
us compare the kinematics when the spectator has a = 1, ps ^ pt and (ps = 180° (in 
plane kinematics) and when a = 1, Ps ^ Pt and (pg = 90° (out of plane kinematics). 



^The reason of this reduction is rather complicated. In the kinematics where = 1 but pt 
is not very small, large effects of rescattering are expected. The deuteron wave function in 
the rescattering amplitude depends on — /cj ~ so the average transferred momenta in the 
rescattering integral are < fcf >~ Pt- The modification due to the A4 and A3 factors leads 
to the condition that a maximal rescattering occurs at a, = 1 + and a. = 1 + Pt— 

for the amplitudes Fb and respectively. Since for 90° cm pp scattering p^t ~ —Pit these 
conditions separate the kinematical regions where rescattering contributions are maximal and 
therefore the maxima of Ff, and Fc are at different a^. As a result the interference terms between 
contributions of different rescattering amplitudes are reduced. 
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Estimating the contribution of FSI amplitude from rescattering diagrams Fig. lb, c it is 
easy to show that the difference between the ps for eq. ([I3| ) is: 

^p2^i^n plane) _ ^2^(out plane) ^ ^2(^2) ^< ^2 > ^ ^2 > ^27) 

Since the deuteron wave function depends strongly on momentum, the difference in the 
argument of the wave function leads to a significant effect in the dependence on azimuthal 
angle. Therefore the FSI interaction for out-of-plane kinematics is larger than for in-plane 
kinematics. This is illustrated in Fig. 5 where the (ps dependence of transparency T is 
calculated at a = 1 for different values of the spectator transverse momenta pt. 

Let us discuss briefly the reliability of the derived picture within the eikonal approx- 
imation. The effects of the order of (~ ^) are included in the argument of the deuteron 
wave function. In perpendicular kinematics the deuteron wave function is sensitive to the 
momentum p ~ 0. Thus ~ effects are amplified by the steep momentum dependence 
of the deuteron wave function at small momenta. In Ref. |^ we demonstrated that the 
off-shell effects are minimal for the perpendicular kinematics since the s dependence of 
the hard pp scattering in this case is minimal. For the kinematics where soft rescatterings 
dominate (at as ~ 1), the calculations are practically insensitive to relativistic effects in 
the deuteron wave function 0. A practical conclusion is that the kinematical restriction 
a ^ 1 and pt < 350 — 400 MeV/c opens a possibility to investigate reliably effects which 
are sensitive to initial and final state interactions of energetic protons with the slow spec- 
tator neutron. Another important question is the validity of the factorization of hard and 
soft scatterings. We demonstrate in Appendix B that the errors due to this factorization 
are minimal for small values of the spectator transverse momenta and > 1. The nu- 
merical estimates of the errors are shown in Fig. 12, which shows that, in the region of 
0.9 < < 1.1 and pt < 400 MeV/c, the factorization approximation is valid within an 
accuracy of better than 15% for energies pi > 6 GeV/c. The accuracy improves at higher 
incident energies (as shown in Fig. 12, for calculations at pi = 6 GeV/c and 15 GeV/c.) 



B. Implication for color transparency 

The aim of the paper is to explore whether the hard exclusive d{p, 2p)n reaction can 
initiate effects of color transparency. The main idea is that at the point of interaction 
the hadrons of the hard elastic scattering are in a "point like" configuration (PLC) 



whose subsequent strong interaction is reduced. p2| , p3| ,|3|,|2T| . As a result soft interactions 
shortly before and after the hard collision will be lower than the usual strong interaction 
of hadrons. Since a PLC, produced in the hard process, is not an eigenstate of the QCD 
Hamiltonian but rather a wave-packet of destructively interfering mass states, it has 
to evolve eventually into a final hadron state. Due to time dilation, the characteristic 
distances for the evolution of a PLC to the normal hadronic state increases with the total 
energy of the PLC. 
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The exclusive d{p, 2p)n reaction has several features which are very sensitive to CT. 
Since ISI/FSI occur at internucleon distances in the deuteron of about ~ 1 Fm p|, one 
can use the spectator neutron to tag the PLC at an early stage of evolution thus reducing 
expansion effects. As a result, large CT effects can be expected even at comparatively 
low energies. Large spectator momenta Ps with ~ 1 ensure small impact parameters 
for soft pn rescatterings. This will enhance the ISI/FSI contribution to the d{p, 2p)n 
cross section. 

We calculate T of eq.(0) in a model which accounts for PLC formation and its time 
development in the framework of the quantum diffusion model (QDM) We note that 
CT in the resonance basis representation of PLC |Tl[| predicts rather similar effects 0. 



We follow the procedure described in ref. p5| to calculate the pn scattering amplitude 



and profile function of eq.(|TT])) within QDM The pn scattering amplitude f^^ in 
eq.(|ll]), (p^ , (|15D , and ( pTsD is replaced by a position dependent one: 

where GN(t) is the Sachs form factor, atot{z,t) is the effective total cross section of the 



interaction of the PLC at distance z from the interaction point. According to ref. |llO| : 



where Ih = 2p/AM^, with AM^ = 0.7 — l.lGeV^. Here {rt{t)'^) represents the transverse 
size of the initially produced configuration. Theoretical analysis of realistic models of a 
nucleon indicate that this size is rather small even for \t\ > 1.5 GeV^. Any 



effects of interplay between large and small size configurations can be included in a 
rather straightforward manner P,^. However our aim here is primarily to demonstrate 
sensitivity to CT. 

The QDM calculation should be compared with the elastic eikonal calculation of the 
previous section, where pn soft rescattering was taken to be practically energy indepen- 
dent (see eq.(^Bp). In order to emphasize any color transparency effects, we calculate 
again the dependence of T on the spectator transverse momentum, azimuthal angle and 
incoming proton momentum, since we showed in Sec. 3.1 that all these quantities are 
rather sensitive to the soft pn interactions. In Fig.6, we present the pt dependence of T 
(eq.(^)) calculated in the eikonal model with eq. (p6D and in the QDM for the rescattering 



amplitude with eqs.(^,^). The shaded areas correspond to the range of the parameter 
AM^ = 0.7 — 1.1 GeV^ which characterizes the scale of excitation energies in PLC and 
which controls the distance over which the PLC evolves to the normal hadronic state. 
Large values of AM^ correspond to small deviations from the elastic eikonal prediction. 
Fig.6 also shows that, at fixed initial energy, CT effects become more prominent with in- 
creasing Pt. This can be understood by the fact that the contribution of the lA amplitude 
is small and the soft rescatterings occur at small internucleon distances leading to some 
suppression in expansion of PLC. Also, higher order rescatterings are more important 
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at larger pf. They are proportional to higher powers of a^°^{z,t) and therefore should 
be more suppressed by the onset of CT. The dependence of T on the azimuthal angle 
of the spectator, in the a = 1 kinematics, is shown in Fig. 7. The ISI (Fig. Id) does not 
contribute to the (ps dependence because we chose the z axis (lab) in the pi direction. We 
see in Fig. 7 that the onset of CT in the (ps dependence requires larger projectile energies. 
This is because the CT impacts the FSI only at energies Ps ~ P4 ~ Pi/2. At energies 
(~ 6 GeV), before the onset of CT, the dependence can be used to check predictions 
of the elastic eikonal approximation. 

To discuss the energy dependence of T we include in our consideration also the models 
which accounts for both the PLC and large (Blob) size configurations in hard scattering 



pp amplitude. We follow to the prescription of Ref. and represent the hard pp 
scattering amplitude as: 

A^-'^(s, t) = A^^^(s, t) + A^^^(s, t). (30) 

A^^^{s,t) is the component of hard scattering amplitude which produces the PLC and 
leads to CT. A^^'-^i^s^t) corresponds to the large size (soft) component of pp scatter- 
ing amplitude and has a cross section comparable to atotai- The source of such a soft 
component could be either the opening of cc channels at energies near threshold 
or the presence of Landshoff processes in hard pp scattering Implications of both 



mechanisms for CT were discussed in details in Ref. ||26|, where PLC expansion effects 
are accounted for within the multi resonance expansion model of CT. In the present 
calculation we describe the expansion effects within QDM using the same Ralston-Pire 
parameterization for A^^'-^{s^t) 0. The numerical results of the multi resonance model 
and QDM for CT in {p,2p), and {e,e'p) reactions are very close 0. In Fig.8 we present 
the energy dependence of T at two values of Pst = 0.2 GeV/c and Pst = 0.4 GeV/c. It 
shows that the higher the spectator transverse momentum, the larger is the CT effect. 
The increase of the incoming energy of the proton diminishes the sensitivity of T to 
AM^, due to reduced sensitivity to the expansion. The model which accounts for the 
interference between large and small size components of hard pp scattering, as expected, 
reveal oscillation with the increase of projectile energy. New feature is that amplitude 
of these oscillations increase with pt and the phases of oscillations are opposite for kine- 
matics dominated by screening {pt = 0.2 GeV/c) and by rescattering {pt = 0.4 GeV/c) 
effects of ISI/FSI. 

As we have seen in Sec. 3.1, T < 1 at pst ^ 0.2GeV/c, and T > 1 at p^t ~ 0.3 GeV/c, 
(see, Figs.3,6,8). Since, by definition, T = 1 means complete transparency (i.e. no 
ISI/FSI) CT produces opposite trends in the regions where T < 1 and T > 1 (see e.g. 
Fig. 6b). This property calls for a more sensitive quantity to characterize CT P,p5|. We 
define the ratio: 

fi=^. (31) 



where p^^^ and p^^^ are such that T{p^^^) > 1 and T{p^^^) < 1. This ratio is more sensitive 
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to CT due to different trends of CT effects in the numerator and denominator of R 
while it will be less sensitive to uncertainties in the theoretical calculations. In Fig. 9 we 
present the energy dependence of R, which shows deviations of QDM prediction from 
the eikonal approximation, in the range of the initial proton momenta > 6 GeV/ c, by a 
factor 2 — 3 due to CT. Because of opposite phases of energy oscillations within the models 
which accounts for the interference between BLC and PLC (Fig.8), R reveals profound 
oscillations in Fig. 9. The dependence of R on the azimuthal angle of the spectator - 
Fig. 10 at fixed energies of initial proton could be used as a complementary method in 
the study of CT. Due to large FSI in the out-of-plane kinematics, the QDM prediction of 
CT at pi = 15 GeV/ c could be enhanced by the factor of 5 — 6 over the eikonal results. 



It was shown in Ref. that lA is strongly suppressed for inelastic recoil final 



states in the deuteron fragmentation region. This was due to the fact that under those 
circumstances the lA is controlled by the inelastic component of the deuteron ground 
state, which is negligible for recoil momenta < 700 MeV/c. Thus the cross section 
is dominated by the rescattering diagrams. Therefore the CT will manifest itself by 
a reduced production of low momentum inelastic states in the deuteron fragmentation 
region at increasing projectile energies. In Fig. 11 we present the calculation of the ratio 
of the cross section of quasielastic d{p, 2p)n reaction to that of d{p, 2p)A^*(1580) reaction. 
The calculations are done in the QDM and the elastic eikonal approximation [^. Fig. 11 
reveals a strong sensitivity of the ratio to CT. For incident momenta from 6 — 18GeV/c 
the eikonal approximation shows a drop in the ratio of 10%, while the inclusion of CT 
leads to a decrease by a factor of 3. This suggests a very simple method for studying CT 
effects i.e. the measurement of the energy dependence of quasielastic and inelastic rates 
in the deuteron fragmentation region. 



IV. SUMMARY 

We calculated the cross section of hard exclusive d{p, 2p)n reaction taking into account 
the initial and final state interactions, within the elastic eikonal approximation. Analysis 
of Feynman diagrams in the relativistic domain produced significant effects beyond the 
conventional nonrelativistic Glauber type formulae. The results can be applied to any 
nucleus. The deuteron calculations were done analytically using a well known form of 
the deuteron wave function. 

The calculation produced a diffractive pattern in the cross section due to ISI/FSI. 
In the kinematics where ~ 1, Pst < 0.2 GeV/c the ISI/FSI suppress the cross section 
with respect to lA, while at ~ 1 Pst > 0.3 GeV/c the ISI/FSI increases the cross 
section. Another result is the prediction of an azimuthal angle dependence due to FSI. 
This is a result of including small but finite angles of the scattered protons into the soft 
rescattering amplitude. We estimated the validity of the factorization approximation 
within the elastic eikonal approach. We found that the factorization approximation is 
valid for ~ 1 to better than 15% for p^t < 0.4 GeV/c. 
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Since ISI/FSI dominate the cross section in d(j), 2p)n, the CT phenomenon will be best 
observed in this reaction. We calculated the transparency T with the quantum diffusion 
model and model where the interference between small and large size components of large 
angle pp scattering are taken into account explicitly. Exploiting the fact that CT leads to 
opposite effects in the kinematical regions where T < 1 and T > 1, we introduced the ratio 
R = ^1^^, which reveals extra sensitivity to CT effect. At presently accessible energies, 
the predictions of energy and angular dependences of R differ by as much as a factor of 
2-6 between the eikonal and CT calculations. The model with the interference of large 
and small size configurations predicts profound oscillations for the energy dependence of 



An additional possibility to search for CT was found. It consists of measuring the 
energy dependence of the ratio of the production rates of elastic and inelastic recoil states 
in the deuteron fragmentation region. While conventional calculation predict practically 
no energy dependence, the CT predicts a significant drop of this ratio with increasing 
the projectile energies. 
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APPENDIX A: ANALYTIC CALCULATION OR RESCATTERING 

AMPLITUDE 

We calculate the rescattering amplitude in eq.(|^) by the method described in ref. 
using the deuteron wave function in momentum space, defined as [|T^: 



R. 



V. ACKNOWLEDGMENTS 




(Al) 



where is the deuteron spin function and 
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SiPz,Pt) 



(A2) 



where Up, Up Pauli matrices. The functions u{p) and w{p) are the radial wave functions 
of S- and D- states, respectively and they can be written as ||T2|,p!3|: 



9 I 2 ' 

+ rrij 



w 



(p) = E 



p"^ + W?j 



(A3) 



where I] = dj = 0, which guarantees that u{p),w{p) ~ ^ at large p and = 

to provide w{p = 0) = 0. Insertion of eg. (|A3|) into the eq.(|^ gives: 



X 



dp's 



+ 



(27r) \ pf^ + m| _^ ^2 



]S{Psz,Pst) 



P'sz - Psz + A4 + ie 



(A4) 



Substituting + = (p^^ + iyrnj + p'st){p'sz ^ "^j + PS) '^'^^ "^^^ perform the inte- 
gration over p'g^ by closing the contour in the upper p'^^ complex semi-plane. Note that 
the dependence of the tensor function S{j)) will not introduce a new singularity, since 
w{p = 0) = 0. Setting the residue at the point p'^.^ = iJm?, + p'st we obtain: 



9 



2iJpf + V 8 



S{iJp? + mlp'J 



2iJpf + m| 



i\lpf + im?j + A4 



(A5) 



After regrouping of the real and imaginary parts, the above equation can be rewritten 
as: 



4i 



(2^ 



(A6) 



where Ps = PsiPsz — ^iiPst — kt), ip'^ is the wave function defined in eq. (|A1|) and ijj'^ is 
defined as: 

i^"{p) = { ui{p)pz + ^^^^^S{p,,p,) + ^ [S{p.,pt) - S{0,pt)]] x", (AT) 



Vsp. 



where 
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j yPt+fn'jip'^ + m'j) 



T \lpi + rn?j{p'^ + 



Note that the last term in eq. (|A7|) does not have a singularity at = since {S{pz,Pt) 
S{0,pt)) 



APPENDIX B: VALIDITY OF THE FACTORIZATION APPROXIMATION 



In this Appendix, we check the validity of factorization of the hard NN scattering 
amplitude from the integrals over the soft rescattering. In the rescattering amplitudes 
(eqs.(|12D, (0), (|15D, (|0D, (P)) s' and t' which enter in the hard amplitude A'^^''''^{s' ,t' 
differ from the measured s 
we have: 



(ps +P4)^ and t = (pi —psY- For the rescattering diagrams 



s' = s 



t' 
t' 
t' 



(Pi -P3- kf 
(Pi -P3- kf 



(Fig.lb) 
(Fig.lc) 
(Fig.lc) 



(Bl) 



where k is the transferred momentum in soft rescattering. To estimate the difference 
between s', t' and s, t it is important to take into account the specific feature of soft, high 
energy collisions, i.e. that the transferred momenta are predominantly transverse with 
respect to the momentum of fast projectile i.e. pjk = 0,( j = 1, 2, 3). The effects neglected 
within the factorization approximation can be estimated by using the hard pp scattering 
amplitude in the form A^^^s^t) ~ ^0.8^3.2 (c.f. Ref. | p8|| ) . With this parameterization we 
obtain: 



X 



A{s,t) 



1 



4s 



2-as 



1 _ 0.8<fc2> i{'2-as)-^Jas{2-as)-l \ / 
^ 4s 2-as J ^ 



12.8<fc2> 



(Fig.lb) 
iM^) (Fig.lc) 
(Fig. Id) 



fB2) 



In deriving eq.( p2|) we use eq.(|^) for the scattering angles in hard pp scattering and 
average over the direction of transferred momenta at soft rescatterings (as a result the 
terms proportional to k vanished). Using for = 1, the relation < k"^ >^ p1^ in high 
transverse momentum of the spectator nucleon we estimate ~ 15% at pi = 6 GeV/c 

and ~ 6% at pi = 15 GeV/c. The overall effect in the full amplitude is smaller, because 
the contribution of the impulse approximation (diagram of Fig. la) does not contribute 
to the error. Eq.( p2D shows that uncertainties due to the factorization approximation 
are minimal at small transverse spectator momenta (due to the relation < k"^ >~ p^J. 
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Also, for fixed spectator momenta, the uncertainties are smaller at < 1. Despite the 
fact that the ratio in eg. (|B2| ) increases at as > 1 the error due to factorization in the 
overall amplitude will still be small since, in these cases, the lA amplitude dominates 
(see e.g. Fig. 12 below and Ref. [Q]). The reason for the dominance of lA in > 1 
kinematics is that, according to eq.(^, the argument of the deuteron wave function in lA 
is smaller than the one in rescattering amplitude. The error due to factorization in the 
double rescattering diagrams of Fig.lg,h is even less important since, in our kinematics, 
they are just a correction to single rescatterings (see Fig. 2). 

To demonstrate numerically the accuracy of factorization, we use two models to de- 
scribe the APP(s, t) amphtude. In the first model we use the fit to hard exclusive pp cross 
section ~ s~^'^[l — cos"^ {6cm)~'^'^]f {t, s), where f{t, s) is a slowly varying function of s and 
t, (see for details Refs. and 10]),) for the cm. scattering angles 9c.m. > 60°, and we 
assume a dependence for smaller angles. In the second model we used a power 

law dependence ~ s~^°[l — cos'^{9cm)~'^'^]f(t/s) down to \t\ ^ IGeV"^ and for smaller t 
we used the ordinary diffraction formulae, as in eq.(^). The results are given in Fig. 12, 
which shows that in the region of 0.9 < < 1.1 and pt < 400 MeV/c the factorization 
approximation is valid within < 15% at pi = 6 GeV/c and even less for higher energies. 
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FIGURES 




FIG. 1. Feynman diagrams of the eikonal approximation for d{p,pp)n scattering. The 
dashed lines describe the amplitude of NN scattering, the full circles represent the hard pp 
scattering amplitude. 
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Pj = 6 GeV , (^s=180' 
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FIG. 2. Transparency T as a function of angle Og between the momentum of spectator 
neutron and projectile proton. The dashed lines are for the conventional Glauber approximation 
where A3 = A4 = Ai = in the rescattering amplitude. The dash-dotted are the calculations 
which include the relativistic correction {Eg — m)^^ in A3 and A4. The solid lines include 
the transverse momentum of the protons in A3 and A4 (e.g. eq.(6)). The dotted lines are the 
same as the solid lines, but without double rescattering terms. 
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FIG. 3. The pt dependence of the transparency T at ag = 1. The dashed line corresponds to 
Pi = 6 GeV/c and solid line to pi = 15 GeV/c. 
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FIG. 4. T/ie pi dependence of the transparency T at ag = 1 for various 
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FIG. 6. The pt dependence of T at as = 1. The solid line is for the elastic eikonal 
approximation which neglects color transparency effects. The shaded area corresponds to T 
calculated within the quantum diffusion model of CT. Dashed and dash-dotted curves correspond 
to QDM calculation with AM^ = 0.7 and AM^ = 1.1 GeV'^ respectively. 
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FIG. 7. T/ie 4>s dependence of T at pt = 0.2 GeV/c, pt = 0.3 GeV/c and pt = 0.4 Gel^/c 
and as = 1. The solid line is the elastic eikonal approximation. The shaded area presents T 
including CT within QDM. The parameter AM^ lies in the range 0.7 — 1.1 GeV^. 
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FIG. 8. T as a function of projectile momentum pi, atpt = 0.2 GeV/c andpt = 0.4 GeV/c. 
The solid line is the elastic eikonal approximation. The shaded area presents T including QDM 
prediction of CT within quantum, diffusion model with the parameter AM^ within the range 
0.7 - 1.1 GeV^. Dash-doted line - QDM prediction with AM^ = 0.7 GeV^, dotted line - 
AM^ = 1.1 GeV^. Dashed line corresponds to prediction CT within the model which accounts 
for interference between large and small size configurations in pp hard scattering amplitude. 
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FIG. 9. The pi dependence of the ratio R = at as = 1. The solid line is the elastic 

eikonal approximation. The shaded area is the the R calculated within quantum diffusion model 
of CT with AM^ parameter in the range 0.7-1.1 Gey2. The definitions of curves are same 
as in Fig. 8 
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FIG. 10. The ratio R = at as = 1 as a function of the spectator azimuthal angle 

(f)s The solid line is the elastic eikonal approximation. The shaded area includes CT, calculated 
in the QDM, with AM^ parameter in the range 0.7 — 1.1 GeV^ . 
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FIG. 11. The energy dependence of the ratio of the cross sections for the reactions d{p, 2p)n 
and d{p,2p)N* . The calculation was done for the elastic eikonal approximation (dashed line) 
and within the QDM model of CT (solid line), with AM^ = 0.7 GeV^. 
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FIG. 12. T as a function of ag for different pt- The solid lines is for the complete 
factorization, the dashed and dash-dotted lines are the two calculations which do not have the 
factorization as explained in Appendix B. 
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